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Abstract 

The Lyapunov stability is established for the iV-soliton solutions in the Lax hierarchy 
of the Benjamin-Ono (BO) equation. We characterize the iV-soliton profiles as critical 
points of certain Lyapunov functional. By using several results derived by the inverse 
scattering transform of the BO equation, we demonstrate the convexity of the Lyapunov 
functional when evaluated at the iV-soliton profiles. From this fact, we deduce that the 

^ 1 N-soliton solutions are energetically stable. 
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Hu{x,t) = -P [°° ^lUdy, (Life) 

7T J-oo V — X 



I. INTRODUCTION 

The Benjamin-Ono (BO) equation describes the unidirectional propagation of long 
internal waves in stratified fluids of great depth. It may be written in an appropriate 
dimensionless form as 

u t + 2uu x + Hu xx = 0. (1.1a) 
Here, u = u(x,t) represents the amplitude of wave, H is the Hilbert transform given by 

u{y,t) 

-oo y — x 

and the subscripts t and x appended to u denote partial differentiation. The BO equation 
can be written as an infinite-dimensional completely integrable Hamiltonian dynamical 
system. 1 ' 2 A common feature of integrable evolution equations is the existence of an infinite 
sequence of conservation laws. The Lax hierarchy of the BO equation is generated by the 
conservation laws, which we shall denote by I n (n = 2,3,4...). The first three of I n read 

h = o / u2dx i ( L2a ) 

2 J— oo 

/°° / 1 1 \ 

oo b" 3 + 2 uHu v dx ' (L26) 

3 1 \ 
u A + -u 2 Hu x + -u\ I dx. (1.2c) 

4 2 / 

In (1.2), the mass conservation has been excluded since it is irrelevant in the present 
analysis. The BO hierarchy is defined by the following nonlinear evolution equations 

where 5/8u is the variational derivative defined by 

d f°° 51 

fc In+2{u + €v) ^ = L^(xy 

When n — 1, (1.3) becomes the BO equation (1.1) with the identification t\ — t while 
when n — 2, it yields the first higher-order BO equation 3 

3 „ 3 



/oo (5/ n _(_2 
j^r--v(x)dx. (1.4) 



( 3 3 \ 

u t2 = [u 3 + -uHu x + -H(uu x ) - u xx j . (1.5) 



Note that the first member of (1.3) reduces simply to a linear equation u to — u x . As 
will be shown below, all the members of the BO hierarchy exhibit the iV-soliton solution 
characterized by the 2N parameters a,j and Xjo(j = 1,2, ...,N) where N is an arbitrary 
positive integer: 

u — un(x — xi, x — X2, x — xn). (1.6a) 



Here 
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Xj = E(- 1 )* +1 ^- a i*- + ^o, (j = 1, 2, .., N), (1.66) 

dj are amplitude parameters satisfying the conditions aj > 0, aj ^ for j ^ k(j, k = 
1, 2, N) and Xj are arbitrary phase parameters. Explicitly, Ujy has a simple expression 
in terms of a tau function / 

d f 

un = «7^hi — , / = det F, (1.7a) 
where F = (fjk)i<j,k<N is an N x N matrix with elements 

fjk = \ x — Xj -\ )5 jk (1 - 5 jk ). (1.76) 

V a j ) a j ~~ a k 

Here, /* is a complex conjugate of / and 5jk is Kronecker's delta. In particular, for N = 1, 
(1.7) represents the 1-soliton solution with a Lorenzian profile 

2ai , . 

Ul = aUx- Xl y + i- (L8) 

A direct proof of (1.7) using an elementary theory of determinants will be presented in 
Appendix A.@ 

The definition of the stability of solitons may be classified according to the follow- 
ing three categories: i) linear (or spectral) stability, ii) energetic stability, iii) nonlinear 
stability. The energetic stability implies that the second variation of certain Lyapunov 
functional becomes strictly positive when evaluated at the soliton solutions. It would 
also lead to the linear stability since the second variation is preserved for the linearized 
equation. In order to extend the energetic stability to the nonlinear stability which deals 
with small but finite amplitude perturbations, one must take into account higher-order 
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nonlinear terms neglected in evaluating the Lyapunov functional and this makes the anal- 
ysis more difficult. In accordance with the above classification of the stability, we shall 
briefly review some known results associated with the stability characteristics of the BO 
solitons. The linear stability of the BO 1-soliton solution has been proved by solving the 
eigenvalue problem associated with the linearized BO equation. 4 A subsequent nonlinear 
analysis shows that the soliton is also stable against small but finite perturbations. 5 As for 
the general iV-soliton solution, its linear stability characteristic has been established by 
solving explicitly the initial value problem of the linearlized BO equation and investigating 
the large-time asymptotic of the solution. 6 ' 7 In the process, the completeness relation for 
the eigenfunctions of the BO equation linearized around the iV-soliton solution has played 
a central role. The recent study demonstrates the orbital stability of the 2-soliton solution 
in which the stability problem has been settled based on the Lyapunov method combined 
with the spectral analysis of the operators associated with the linearized BO equation. 8 
The approach used in this paper originates from the stability analysis of the multisoliton 
solutions of the Korteweg-de Vries (KdV) equation by means of the constrained varia- 
tional principle. 9 See also an analogous work dealing with the spectral stability of the 
multisoliton solutions in the KdV hierarchy. 10 All the works mentioned above are con- 
cerned with the stability of solitons for the BO equation. The stability characteristics of 
solitons in the BO hierarchy have not been considered as yet. 

The purpose of this paper is to establish the Lyapunov stability of the general N- 
soliton solution (1.6). To be more specific, let us consider the following higher-order BO 
equation which consists of the commuting flows of the BO hierarchy 

d SH N 

where the Lyapunov functional Hn is given by 

N 

H N (u) = I N+2 + Vnln+l, (1.9b) 
n=l 

and \i n are Lagrange multipliers which will be expressed in terms of the elementary sym- 
metric functions of a± : a 2 , ajv- See Sec. Ill for the detail. We define the profile (or 
shape) of the iV-soliton solution by Un = Un(x) = •u JV |t =t 1 =...=o- We observe from (1.7) 
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that the iV-soliton profile has the same functional form for all the members of the hi- 
erarchy, the only difference being the velocities of the solitons. We show that Un is a 
stationary solution of (1.9) decaying at infinity. Namely, Un is realized as a critical point 
of the functional H N . Using (1.9b), this condition can be written as the Euler-Lagrange 
equation 

SIn+2 m-^ SIn+i n , T1 ri in\ 

5u + 2^ l 2n ~^~ = °' at u = Un - ( 1>10 ) 

The Lyapunov stability of Un may characterize Un as a minimal point of the functional 
In+2 subjected to iV constraints 

I n+1 (u) = d n ,(n= 1,2,.. .,N), (1.11) 

where d n are real constants and consequently the second variation of Hn is strictly positive 
at Un- This means that H N is convex at Un, so that the following inequality holds 

H N {U N + ev) - H N {U N ) > 0, (1.12) 

where ev is a perturbation imposed on Un which belongs to certain function space specified 
by the 2N integral conditions. We assume that the L 2 norm of v is finite. The small 
parameter e has been introduced to measure the magnitude of the perturbation. The 
inequality (1.12) shows that the iV-soliton solutions are energetically stable. We give 
a direct proof of (1.12) with the aid of the results obtained by the inverse scattering 
transform (1ST) for the BO equation. 1,2,11,12 In Sec. II, we summarize the background 
results arising from the perturbation theory and the Hamiltonian formulation of the BO 
equation which provide the necessary machinery in carrying out the stability analysis. 
In Sec. Ill, we prove the inequality and hence establish the Lyapunov stability of the 
iV-soliton solutions in the Lax hierarchy of the BO equation. In Appendix A, we present 
a direct proof of the iV-soliton solution (1.7). In Appendix B, we evaluate the number of 
positive eigenvalues of the Hessian matrix associated with H N - 

II. BACKGROUND RESULTS OF 1ST 

The 1ST has been applied successfully to solve the initial value problem of the BO 
equation. 11,12 Furthermore, for real generic potentials it has been used to prove the com- 
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plete integrability of the BO equation. 1 ' 2 Here, we summarize some background results of 
the 1ST necessary for the stability analysis. 
A. Eigenvalue problem 

The eigenvalue problem associated with the 1ST of the BO equation may take the 
form 

#+ + a(0+ - r) = (2.1) 

where + (0~) is the boundary value of the analytic function in the upper(lower)-half 
complex x plane, u is a real potential rapidly decreasing at infinity and A is the eigenvalue 
(or the spectral parameter). We define the two Jost solutions of (2.1) specified by the 
boundary condition as x — > +00 

N(x, A) -> e iXx , N(x, A) -> 1, (2.2) 
and the analogous ones as x — > —00 

M(x, A) -> 1, M(x, A) -> e iAx . (2.3) 

These solutions satisfy the linear integral equations 

N x -i\N = iP+(uN), (2.4a) 

jV x - iAiV = iP+(wiV) - iA, (2.46) 

Mr - iAM = iP+(uM) - iA, (2.4c) 

M x - iXM = iP + (uM) , (2 Ad) 

where P + is the projection operator defined by P + = |(1 — iif). The solutions of (2.4) 
subjected to the boundary conditions (2.2) and (2.3) exist for A > 0. The Jost functions 
M, N and N are then related by 

M = N + (3N, (2.5) 

where (3 is a reflection coefficient. For pure soliton potentials, this reflection coefficient 
vanishes identically. 
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There exists a set of solutions ®j(x) for negative A = Xj(j = 1,2, ...,N) which satisfy 
the equation 

$ jtX - iXjQj = iP+iuQj), (j — 1, 2, N), (2.6a) 
with the boundary conditions 

-> -, x -> +oo, (j = 1, 2, iV). (2.66) 

B. Conservation laws 

It follows from (1.1), (2.4b) and the time evolution equation for N 

N t - 2XN X - iN xx - 2{P + u x )N = 0, (2.7) 

that the quantity j™ 00 u{x,t)N{x,t)dx is conserved in time. Expanding N in inverse 
powers of A 

N=J2 [ \n\ N 1 = l, (2.8a) 

n=0 A 

and substituting (2.8a) into (2.4b), we obtain the following recursion relation that deter- 
mines N n : 

N n+1 = tN n , x + P + (uN n ), n > 1. (2.86) 
The nth conservation law may be taken as 

/oo _ 
uN n dx, (2.9) 
-oo 

where a factor (—1)™ is multiplied for convenience. The first three of /„ except h are 
already given by (1.2). In terms of the scattering data (5 and Xj, I n can be evaluated as 

N 



/„ = (-!)" hn^i-^T- 1 + [ ~^ f ^- 2 P*WmdX , (n = l,2,...). (2.10) 



J'=l 

The first term on the right-hand side of (2.10) is the contribution from solitons and the 
second term comes from radiations. It is important that both contributions are additive. 
A remarkable feature of the conservation laws is that they are in involution, namely 
I n (n = 1,2,...) commute each other in an appropriate Poisson bracket. In particular 
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C. Variational derivatives 

The variational derivatives of the scattering date with respect to the potential are 
calculated explicitly. In developing the Lyapunov stability, we need the formulas of the 
variational derivatives evaluated for the iV-soliton potential u = Un- In particular, the 
following formula plays an important role in our analysis 

(m)^ v -^'' (x)Hx) ' ( ^ 1 ' 2 '---' N) - (2 ' 12) 

Here, the eigenfunction $j corresponding to the discrete spectrum Xj satisfies the system 
of linear algebraic equations 

jv l 

( x _ 7 .)$. + i £ r $fc = 1 ' C? = l,2,...,iV), (2.13) 

fe=i A 3 ~ A k 

where 7, = Xjo + i/(2Xj) and Xj are real constants. Recall that Xj are related to the 
amplitude parameters aj introduced in (1.7) by the relations Xj = —a,j/2(j = 1,2, ...,N). 
Taking account of the fact that the reflection coefficient (3 becomes zero for u — Un, we 
can derive from (2.10) and (2.12) the formula 

=(-l)>-l)f;(-A J -) B - s ^(x)$ i (x), (n = 2,3,...,iV). (2.14) 

\ 0U \ X )J U =U N j=l 

In terms of $j, Un has the following two alternative expressions: 

N 

U N = tJ2(^-®*), (2.15) 
i=i 

jv 1 
3=1 A J 

The positive definiteness of C/jv is obvious from (2.16) since all Aj are negative quantities. 
One can derive (2.16) by using (2.13) and (2.15). The formula (2.16) also follows from 
(1.2a), (2.10) and (2.14). In Appendix A, we show that Un can be rewritten in a compact 
form in terms of a determinant. 

The following relation concerning the variational derivative of j3 with respect to u is 
useful in evaluating the contribution of the continuous part to the functional Hn- 

5 ^ = iM(x,X)N*(x,X). (2.17) 



For the iV-soliton potential u = U N) M reduces to N by (2.5) and (3 = 0. The function 
MN* satisfies the orthogonality conditions 

£^ M(x, X)N*(x, A)^ (^*(x)^(x)) dx = 0, (j = 1,2,..., N). (2.18) 

Finally, we empasize that all the results presented here are obtained through the analysis 
of the spatial part (2.1) of the Lax pair for the BO equation. 

III. LYAPUNOV STABILITY 

A. Variational characterization of the iV-soliton profile 

We first show that the stationary solution Un of the higher-order BO equation (1.9) 
satisfies (1.10) if one prescribes the Lagrange multipliers /i n appropriately. This provides 
a variational characterization of Un- Let tyj = and bj = —\j = %/2. With this 

notation, (1.10) and (2.14) give a linear relation among \F,- 

N N N 

(N+^^bf-^j + £(-l) N -" + Vn = °- 

j=l n=l j=l 

In view of the fact that are functionally independent squared eigenfunctions, 1 ' 2 /i n 
must satisfy the following system of linear algebraic equations: 

N 

£(-l)"-^-Vn = (N + l)bf, (j — 1, 2, N). (3.2) 

n=l 

To solve (3.2), we introduce an N x N matrix V 

V = (v jk )i<j,k<N, v jk = b k ~ l , (3.3) 

and the cofactor of Vjk by 

d\V\ 

V jk = ^LJ, V = det V, (3.4) 
dv jk 

where \V\ is the Vandermonde determinant. Notably, since |V| = Yii<j<k<N(bk — bj) and 
bj 7^ b k for j 7^ k, \V\ never vanishes. This fact will be used essentially in the following 
calculation. It is also convenient to define the polynomials g(x) and g^{x) by 

N N 

g{x) = X[{x-bj) = Y.{-l) s °sX N -\ (3.5) 

j=l s=l 

9 



N N-l 



9k(x) = I] (* - b j) = E (-l) s Vk, s x N -\ (3.6) 



3=1 S=l 



where cr = 1 and a s (l < s < N) are elementary symmetric functions of b±, b 2 , &/v 

AT AT AT 

o"i = E & i> °"2 = E bjb k ,..., a N = Y[bj, (3.7a) 



7 = 1 3,k = l 1 = 1 

(J<fc) 



and <T fc s are given by the relation 

^M = E^(-^r j . (3-76) 
i=o 

Obviously, all Oj are positive quantities since bj > 0(j = 1,2, N). Now, applying 

Cramer's rule to (3.2) with use of the fact \V\ 7^ 0, we find that \i n are determined 
uniquely as 

^ = (-1) JV -" 7V+1EL ^ , (n = 1, 2, AT). (3.8) 
n I V I 

Substituting the formulas 13 

^= ( " 1) ^*f~ w|V| , (fc,n = l,2,...,N), (3.9) 

E^7TT- = ^-^i. (n=l,2,...,JV), (3.10) 
fc= i S'fclOfcJ 

into (3.8), we arrive at a simple expression of ji n 

N + 1 . . . 

A*n = a N - n+1 ,(n = 1,2,...,N). (3.11) 

77/ 

If we use the relations bj = a J /2(j = 1, 2, AT), we can see that \x n are expressed in terms 
of elementary symmetric functions of a±, a 2 , a^. 
B. Stability 

Let us now prove the inequality (1.12) which assures that the functional H N is convex 
at the TV-soliton profile U N . The method used here is based on the ideas developed in a 
recent work 10 on the spectral stability of the iV-soliton solution of the KdV hierarchy as 
well as an earlier work 14 on the algebraic structure of the BO iV-soliton solution. We first 
rewrite (2.10) as 

I n+1 (u) = (-lf +1 j 2 ^X>? + ("l) n+1 r„| , (3.12a) 

10 



where we have put bj = —Xj and 

1 r°° 

r„ = — / A^ 1 /3*(A)/3(A)rfA. (3.126) 

Z7T JO 

Let AQ be the increment of any functional Q(u) around u = Un , i.e., 

AQ = Q(U N + ev) - Q{U N ). (3.13) 

It then follows from the constraints (1.11) that 

A/ n+1 = 0,(n = l,2,..,iV). (3.14) 

We then use (3.12) to rewrite (3.14) in the form 

N 

2nn ]T b^Abj + (-l) n+1 Ar n = 0, (n = 1, 2, N), (3.15) 
i=i 

where we have neglected the higher-order terms (Abj) s (s = 2,3, ...,7V). These relations 
indicate that the increments of soliton amplitudes are balanced with the increments of 
radiations. We recall that f3 = for u — U N and consequently A(f3*f3) = A(3*A(3. This 
leads to the estimates Ar n ~ 0(e 2 ) and Ar n > for all n. The case Ar n = calls a special 
attention and it will considered in detail later. Hence, (3.15) can be solved consistently 
in Abj only if Abj ~ 0(e 2 ). Since by the definition (1.4) 

'5bi 

/ u=Un 

one must impose the integral conditions on the perturbation v(x) 

• x i '5bi 

/ u=Un 

in accordance with the above estimate for Abj. We can see from (2.12), (2.14) together 
with the relations bj = —\j(j = 1,2, ...,N) and \V\ ^ that (3.17) are equivalent to 



Abj = e J°° ( S -^\ v(x)dx + 0(e 2 ), (3.16) 



/_^(ff) v(x)dx = 0, (j = l,2,...,iV), (3.17) 



I- f^r) v(x)dx = 0, (n = 1,2, ...,7V). (3.18) 
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Owing to (3.14), however, these conditions are satisfied automatically. The above obser- 
vations allow us to solve (3.15). Indeed, the solutions are written, with use of Cramer's 
rule, as 

i t n V Ar 

Ab J = ^ |V| ^ ^ (J' = 1.2,...,JV). (3.19) 

It now follows from (1.9b), (3.12) and (3.14) that 

AH N = (-1) N ^(N + VpbfAbj + (-1)^+^+1 1 • (3.20) 

If we substitute (3.19) into (3.20) and use the fomulas (3.9) and (3.10), AH N simplifies 
to 

N 

AH N = (N + 1) ^ -Z=m&r n + Ar N+1 . (3.21) 

n=l n 

Since <7jv_ n+ i > for n = 1,2, ...,N by the definition (3.7), we find that if at least one 
of Ar n is not zero, then AH^ > 0. On the other hand, if all Ar n become zero, then 
AH N = 0. In the latter case, we see from (3.15) that Abj = for all j. This situation 
will happen when the perturbation ev represents the small variation of Un with respect 
to the phase parameters Xjo- Specifically 

^) = £f^o, (3.22) 

j=l OX j0 

where Sxjq are small perturbations of order e. If we impose the following iV integral 
conditions on v(x) in addition to (3.28) 

\ / u=Un 

then the perturbation of the form (3.22) ceases to be admissible, as we shall now demon- 
strate. We first notice that the right-hand side of (3.22) can be expressed in terms of 
the x derivative of (5I n+ i/5u) u= u N . Indeed, we take u = ujy in (1.3) and then put 
to — h — ... = to obtain 

(- 1 '""^rg=|(^) ,.<»=1.2 ^ (3-24) 

j=l J v \ / m=(7jv 
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where we have used (1.6), the definition of U N and bj = aj/2. Thanks to the fact |V| 7^ 0, 
the relations (3.24) can be inverted to give 

An alternative expression of (3.22) follows immediately upon introducing (3.25) into 
(3.22), which reads 

j=l n =l \ \ \ / u =Un 

We observe that this perturbation satisfies the conditions (3.18) by virtue of (2.11). It is 
important that the N x N matrix C = (cjk)i<j,k<N with elments 

Cit= rU^l) U^l) dx, ,3.27) 
7_oo dx V 5u J u=Un dx V 5u J u=Un 

is positive definite and hence |C| 7^ 0. In view of this fact, we deduce from (3.23) and 
(3.26) that 5xj = 0(j = 1,2, ...,N) and consequently v — 0, which implies the assertion 
mentioned above. An additional relation which deserves remark is 



f 

J — < 



-f] v(x)dx = 0, (3.28) 

MJu=U N 



which follows from (2.14), (2.17), (2.18) and (3.26). This leads to the estimates Af3 ~ 
0(e 2 ) and Ar n ~ 0(e 4 ). As a result, the pertubation (3.26) gives rise to higher-order 
contributions to (3.21) which also means that the second variation of Hn turns out to 
be zero. In conclusion, the inequality AH N > holds under the simultaneous conditions 
(3.18) and (3.25) imposed on v(x), which completes the proof of (1.12). The convexity 
of H N implies that the second variation of H N is strictly positive and consequently the 
iV-soliton solutions are energetically stable. 

C. Remark 

In this paper, the convexity of H N has been proved by invoking some results obtained 
by the 1ST of the BO equation. There exists, however another method to establish the 
same convex property without recourse to the 1ST. To illustrate this, we put u(x, t) = 
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Un(x) + ev (x)e xt and linearize (1.9) around U^. The resulting eigenvalue equation can 
be written as 

-^C N v = Xv, (3.29) 
where Cn is a self-adjoint operator. This operator may be defined through the relation 

S 2 H N = — v(x)jCnv(x)cIxi (3.30) 

2, J —CO 

where 5 2 Hn denotes the second variation of H^. Let ti(Cn) be the number of negative 
eigenvalues of Cn and p(Hn) be the number of positive eigenvalues of the Hessian matrix 
defined by 

Hn = (hjk)i<j,k<N, hjk = t: — n — • (3.31) 

Then, under the conditions (3.18) and (3.23) the positivity of 5 2 H N is satified if and only 
if u{Cn) = p(Hn)- The above criterion of the positivity property has been proved in 
Ref. 15 and has been applied to the Lyapunov stability of the A-soliton solution of the 
KdV equation. 9 In particular, the spectral property of the 2Ath-order differential oper- 
ator associated withe the linearized KdV equation has been investigated by extending 
the classical Sturmian theory. See also a related work dealing with the stability of the 
iV-soliton solutions in the KdV hierarchy. 10 In the case of the BO equation, however the 
eigenvalue equation (3.29) is not purely differential equation but actually integrodiffer- 
ential equation since it includes the Hilbert transform. This makes the spectral analysis 
more difficult. Quite recentry, a new method was developed to characterize the spectral 
property of Cn for N = 2. 8 The extension to the general A-soliton solutions of the BO 
equation and its hierarchy is still to be resolved. It is noteworthy that p(Hn) can be eval- 
uated explicitly for the iV-soliton profile U^. This calculation is presented in Appendix 
B. The stability analysis developed in this paper would suggest that u(Cn) is equal to 
p(Hn)- This interesting issue will be persued in a future study. 

APPENDIX A: PROOF OF THE A^-SOLITON SOLUTION 

In this Appendix, we provide a direct proof of the A-soliton solution (1.7) of the nth 
higher-order BO equation (1.3) by means of an elementary theory of determinants. For 
convenience, we write down some basic formulas for determinants upon which our proof 
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relies. Let F be an iV x iV matrix with elements f jk given by (1.7b) and F jk be the 
cofactor of fj k The expansion of \F\ by elements and their cofactors is given by the two 



ways: 

TV 



Y,fikF lk = 5 3l \Fl (Ala) 



k=i 

TV 



J2fjkFji = 5 kl \F\. (Alb) 

3=1 



The folowing formula is a consequence of (Al) 



TV TV 

E Ui + 9k)f jk F jk = + 9j )\F\. (A2) 

j,k=l j=l 



The differential rule applied to the determinant \F\ gives 

TV 
3=1 

I^L = (-!)" E^i, (^36) 

where Cj = (n + l)a™/2 n . To carry out the proof, it is necessary to assign the time 
dependence of the eigenfunction $j for the discrete spectrum Xj. This can be accomplished 
simply by replacing the phase factor jj introduced in (2.13) by Xj which is defined in (1.6). 
We first show that (2.15) can be rewritten in an alternative determinantal form (1.7). The 
solution to (2.13) is found by using Cramer's rule as 

TV jp 
fkj 



*; = < E||t- ^ 



k =i 



We put fj = dj and gj = —dj in (A2) to derive the relation 



It follows from (A3)-(A5) that 



TV TV 

E /> E r ,r (A5) 

3,k=l 3=1 



E^- = i(ln|F|) x . (A6) 

3=1 
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Substituting (A6) and its complex conjugate expression into (2.15), we find that (2.15) 
coincides with (1.7). 

Let us now proceed to the proof of the A-soliton solution. We substitute (1.7) and 
(2.14) into (1.3) and integrate it once with respect to x to recast (1.3) into the form 



N 'a^™- 1 



i 



(1*1*1*1*. - |F||F] t *J/|F|*|F| = (-1)> + 1) £ [-£ ) * (A7) 



where we have used the relation \j = — %/2. The following identity has been established 
by using Jacobi's formula for determinants: 

jwf* _ mu 2 Eww, 0;t=1>2 ,...,„ ) . ^ 



\ Clk a>j j OjjOjk 

Indeed, (A8) coinsides with (A20) in Ref. 7 with the identification / = |F|*,Aj^ = 
F* k ,ipj = $*. If we multiply (A8) with j = k by Cj and sum up with respect to j, we 
obtain 

% N n+ 1 N 

^AW\*Fn ~ W\Fh) = ^rY.*T 1 ^\F\*\F\. (A9) 

A 3=1 / 3=1 

The left-hand side of (A9) is modified further by introducing the formula (A3b) and its 
complex conjugate expression. It leads, after dividing the resultant expressin by |F|*|F|, 
to (A7) and thus completing the proof. 

APPENDIX B: POSITIVE EIGENVALUES OF THE HESSIAN MATRIX H N 

The Hessian matrix is defined by (3.31). It is a real symmetric matrix whose 
elements are calculated explicitly for the A-soliton solution. Indeed, by taking (3 = in 
(2.10), the nth conservation law corresponding to u = un reduces to 

N 

I n = 2v{-l) n Y.Vr\ (6i A,). (Bl) 

i=i 

If we regard H N as a function of Uj(j = 1, 2, N), we obtain from (1.9b) and (1.10) 

dH N 
dm 

Hence 



/ j+1 ,y=l,2,...,JV). (B2) 
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Let P = (pjk)i<j,k<N and Q = (qjk)i<j,k<N be N x N matrices with elements 

p jk = 2n(-iy +1 M~\ (54a) 

_ dfij N + 1 da N - j+1 

Qj k - - : o7 > (#4o) 

db k j db k 

respectively. Note that the right-hand side of (B4b) follows from (3.11). Using the above 
definition, we can rewrite (B3) in the form 

H N = PQ~\ (B5) 

if Q^ 1 exists. To show the nonsingular nature of Q, we use the definition (3.7a) of a^-j+i 
and (B4b) to evaluate the determinant of Q. A simple calculation immediately leads to 

(N + 1) N 

\Q\= { l } ' II (B6) 

l<j<k<N 

Since bj ^ bk for j ^ k, we confirm that \Q\ ^ 0, implying that Q is invertible. 
It now follows from (B5) that 

Q T H N Q = Q T P. (B7) 

In accordance with Sylvester's law of inertia, one can see from (B7) that the number of 
positive eigenvalues of 7ijv coincides with that of Q T P. The latter can be counted easily, 
as we shall now demonstrate. Using (B4), the (j, k) element of Q T P becomes 

N Q 

(Q T P) 3k = 2,(N + 1) E(-l)' +1 ^ ±i ^ 1 - (B8) 

i=i 00 i 

We differentiate (3.5) by bj and then put x = b k to derive the relation 

N r)rr 

E(-l) m %^^ = n (6« - b k ). (B9) 
i=i O0 j i=i 

Introducing (B9) into (B8), we find that Q T P is a diagonal matrix. We can order the 
magnitude of bj as &i > b 2 > ... > b N > without loss of generality. Then, (B8) and 
(B9) indicate that the number of positive eigenvalues of Q T P is equal to where [x] 

denotes the integer part of x. If we take account of (B7) and Sylvester's law of inertia, 
we conclude that p[Hn] 



N+l 
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